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To explain the cosmic speed up, brought to light by the recent SNIa and CMB ob- 
servations, we propose the following: a) In a spacetime endowed with a FRW metric, 
we choose an empirical scale factor that best explains the observations, b) We assume 
a modified gravity, generated by an unspecified field lagrangian, f(R). c) We use the 
adopted empirical scale factor to work back retroactively to obtain f(R), hence the term 
'Inverse /(-R)'. d) Next we consider the classic GR and a conventional FRW universe 
i that, in addition to its known baryonic content, possesses a hypothetical 'Dark Energy' 

(~^) ' component. We compare the two scenarios, and find the density, the pressure, and the 

' equation of the state of the Dark Energy required to make up for the differences between 

f"*"^ , the conventional and the modified GR models. 
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As cosmological standard candles, supernovae type la (SNIa) appear dimmer 
than what one expects from a Cold Dark Matter (CDM) model of the universe 
| 1 | 2 | 3 |_ Thjg observation and other evidences from the Cosmic Microwave Background 
(CMB) m easure ments indicate that the universe is in an acceleration phase of its 
expansion I 4 | 5 | 6 | A conventional CDM scenario does not explain this speed up. 
Some authors have stipulated a dark energy component to make up for whatever 
dynamical effect s that the known energy momentum co ntent of the model leaves 
unaccounted for I 7 | 8 | 9 | 10 | ll | 12 | 13 | 14 | 15 | 16 ll7ll 8 19 20 2T| p^er s have entertained 
alternatives to Einstein's gravitation | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 |30l Y et a third school 
have resorted to inhomogeneous FRW universes to explain the dilemma Since 
all these approaches attempt to answer the same question, all should be equivalent, 
and there should be a way to translate one language to the other. 

Here, we are concerned with the 'dark energy' and 'alternative gravitation' sce- 
narios. We suggest to begin with a Freidman- Robertson- Walker (FRW) universe, 
to choose its scale factor, a(t), in a way that best explains the available observa- 
tions, and to work out the dynamics of the spacetime. Next, to write down the 
field equations for a modified f(R) gravitation, and knowing a(t), to solve for f(R). 
Finally, to attribute whatever deviations from the conventional FRW results there 
is, to a dark energy field, and to obtain its density, pressure, equation of state, etc. 

The choice of the scale factor 

In a FRW metric, ds 2 — —dt 2 + a(t) 2 dx. 2 , single term scale factors of the form 
a oc t@ lead to constant deceleration parameters, q = —da/d 2 = (1 — /?)//?, and do 
not serve the purpose. We propose the following two-term ansatz 

a(t) = -^—(t/to) 2 / 3 \l + P (tAo) 2Q/3 l , (1) 
1+p L J 

where to is the age of the universe, and a and p are the free parameters of the model, 

to be adjusted to ensure compatibility of the emerging results with observations. 
The factor (1 +p)~ x is introduced to have a(to) = 1. By letting either a or p tend to 
zero, one recovers the standard CDM universe. Hereafter, for economy in writing, 
we will use the time parameter r = p(t/to) 2a ^ 3 instead of the conventional time t. 
From Eq. JT]) one finds 

q = \ [1 - (1 + a)(2a - l)r] [1 + r] [1 + (1 + a)^ 2 (2) 



2 

2 

3t y p' 



H=^-C-)- 3 ' 2a [l + (l + a)r][l + rr\ (3) 



TZ = t 2 R = 6H 2 (l- q ) = -(-)-^ a (4) 

3 p 

x [1 + (2 + 5a + 2a 2 )r + (1 + 5a + 4a 2 )r 2 ] [1 + r] -2 . 

For a > 1/2, q can become negative and remain nonsingular throughout. Transition 
from a decelerated phase of expansion to an accelerated one takes place at r trans = 
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[(a + l)(2a - l)]" 1 or i trans - [(a + l)(2a - l)p]" 3/2a to < t . For -1 < a < 1/2, 
q is always positive and nonsingular. For a < — 1, q can become negative but also 
singular. The last two possibilities are discarded. For all values of a and p, 

limit q(t->0) = i, 

limit q(t -> oo) = (1 - 2a)/2(l + a). (5) 

In Fig- (JTJ) we have plotted q(t) versus the redshift, z = a(to)/a(t) — 1, for several 
values of a and p = 1/3. 




Fig. 1. Plot of q(t) versus the redshift, z = a(t )/a(t) - 1, for a = 0,1, 2, 3,4; and p = 1/3. The 
case a = gives the classic value, q = 2/3. As a increases from 1 to 4, transition to the accelerated 
phase of expansion moves from later to earlier epochs, from smaller z's to larger one. 



As to H and R, both remain positive for all times. Both tend to oo as r — * 
and decease monotonically to as r — > oo. They exhibit a normal behavior in the 
neighborhood of T trans = [(1 + a) (2a — l)] -1 . 

Equation ([3|), written for the present epoch, reveals a constraint on a and p, 
that should be observed in the final design of the model. Thus, 



[l + (l + a)p]/[H-p] = |ff *o«|. 



(6) 
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Fig. 2. The distance modulus of the SNIa Gold's sample versus redshifts, black circles; And the 
plot of Eq. (0 with the scale factor of Eq. JTJ, solid line. Parameters of the model are a = 2 and 
P = 1/3. 

Empirical values of a & p 

The distance modulus (corrected for the reddening) and the (dimensionless) lumi- 
nosity distance, Dl(z; a,p), of supernovae are related as 

(i = m-M = 5\ogD L (z\a,p) + 25, (7) 

D L (z;a,p) = (1 + z) J dzH(z; a,p) _1 . 

In Fig. @, the observed distance modulus of 157 SNIa of Gold's sample are plot- 
ted versus the redshifts. The solid curve is the plot of equation ([7]), in which, in 
compliance with the constraint of Eq. ([5]), we have chosen 

a = 2, p=l/3. (8) 

The fit to the data points is adequate for our purpose, though the parameters can 
be refined to optimize the fit. These numbers give itrans = 0A3t Q and z trans = 1.14. 
On adopting Ho w 70km sec^ 1 Mpc -1 , one obtains an age of to « 12.4Gyr for the 
universe. 

An inverse f(R) way out 

We begin with a modified field equation generated by an, as yet, unspecified field 
lagrangian, f(R), 

R»u - \r9h» = ^p(f~ R p )9^ (9) 
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where F(R) — df(R)/dR, and k = 8ttG. For a universe of FRW type, filled with a 
perfect fluid of density p rn and pressure p rn , Eq. ([9]) and the equation of continuity 
reduce to 

3HF + 3H 2 F- l(RF-f)-K Pm = 0, (10) 

F-HF + 2HF + K {p m + Pm ) = 0, (11) 
Pm + 3H(p m +p m ) = 0, (12) 

We further neglect the pressure and integrate Eq. (TT2"j) to obtain p m (t) — p a(t)~ 3 . 
Next we substitute for p m in Eq. (jlip. assume a = 2, change the time variable to 
t = p(t/t ) 4 / 3 , and find 



(l + r) 3 r 2 F" 



- i (l + r) 2 (l + 5r)rF' 



(13) 



3r 2 )F+- =0, 



where the '/ ' now stands for d/dr, and we have, arbitrarily, put the dimensionless 
constant (3/4) 3 (l +p) 3 to 2 npo, that appear in the course of mathematical manipula- 
tions equal to one. We will shortly discuss the numerical solution of Eq. (| 13|) . Some 
general remarks on its asymptotic behavior, however, are instructive. As r — > we 
find 

23 9 



F(r) = 



7 

-0.44 



G 



(14) 



d r- u - 44 F 1 (r)+ C2 t^P 2 (t), t 



0. 



where Pi and Pi are calculable polynomials in r, and begin with term 1; cx and C2 are 
constants of integration to be obtained from boundary conditions. The exponents, 



-0.44 and 1.7 are approximate solutions of the indicial equation, s 2 



0. As 



t — ► 0, F diverges to oo or converge to on account of one or the other term. This 
feature makes the solutions sensitive to a CDM type boundary conditions of the 
form, -F^Tinitiai) = 1 and -F'(ri n iti a i) = 0. Presently we have no basis, observational 
or otherwise, to make an intelligent guess as to what the appropriate boundary 
conditions are. For the sake of argument, however, we have adopt c\ = c 2 = 0, and 
kept only the proper solution of Eq. (fT2|) . With F(t) known, Eq. (fT0|) becomes an 
algebraic equation to calculate /(r). The numerically calculated solutions of F(t), 
/(t) and 1Z = io 2 i?(r) are plotted in Fig. (J3j> . Elimination of r in favor R provides 
F(R) and f(R) 

It is instructive to examine the asymptotic behavior of F(R) and f(R) analyti- 
cally. In the limit of small and large r's, corresponding to large and small i?'s, one 
finds 



F(R) = 



6 



n 



n 
T 



5 / TZ 
28^3 \p 



2/3 



n -> o. 



(15) 



(16) 
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Fig. 3. F(t), proper solution of Eq. H13I I (xlOO), dot- dashed line; /(t), dashed line; and 1Z(t), 
solid line; a = 2, p = 1/3. 



Note that F(R) is a dimensionless scalar as it should be. With F(R) 
known, it is a matter of simple integration to obtain f(R). Thus, 



f(R) 



K 

*7 



+ ■ 



11 



i 



15 



112^/3 \p 



K - 

2/3' 



n-> o. 



df/dR 

(17) 
(18) 



At early epochs (large -ff's), the spacetime approaches the conventional FRW uni- 
verse with the classical GR while for the later times we have a positive acceleration 
universe. Another point about this action is that for R — > in the solar system 
scales, /(0) = and /'(0) = 0. In this case we will have standard GR equation 
and f(R) evades from the solar system test. Recently this type of models have been 
studied by introducing action in the form of f(R) = R + fi(R), where for R = in 
the solar system, /(0) = and for larger i?'s, in cosmolog ical sc ales and inside the 
large scale structures the action reduces to f(R) = R — A ^ l ^ i 
In the remaining range of 1Z, integration is done numerically and the results are 
plotted in Fig. g]). 



Dark Energy equivalent 

Instead of the modified gravitation considered above, let us assume a classic FRW 
universe. That is, let f(R) = R and F = 1. Let this universe, however, have a 
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40 




Fig. 4. Numerical plot of f(R) versus R. r is eliminated between /(r) and R(t); a = 2, p = 1/3. 

'Dark Energy' component, in addition to its conventional baryonic content. The 
counterparts of Eqs. (TO]) , and (fTT|) will be 

3H 2 - K (Pde + Pm) = 0, (19) 

2H + k (p de + p m ) + k (p de + p m ) = 0. (20) 
Subtracting Eq. (10]) from (19]), and Eq. JTT]) from (20]) gives 

K Pde = 3H 2 {l-F)-3HF+^(RF-f) 1 (21) 
K( Pde + Pde ) = F-HF-2H(l-F), or 

K Pde = F + 2HF -H 2 (l- 2q)(l - F) 

-\{RF-f). (22) 

The equation of state for the dark energy is obtained by eliminating t, implicit in 
F and H, between Eqs. ([2~T]) and ([22]) . This is done numerically and w = p de /pde 
as a function of the redshift is plotted in Fig. (5]). 

Concluding remarks 

The algorithm of Fig. [5] summarizes the path we have followed in this letter. We 
have resorted to the SNIa observations to design an empirical FRW metric that 
allows the model universe to transit from a phase of decelerated expansion at early 
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Fig. 5. Equation of State of the Dark Energy: Numerical plot of w = Pde/pde versus the redshift 

epochs to an accelerated one at later times. The spacetime is almost a CDM model 
and the gravity is almost the classic GR at very early times, but evolves away in 
course of time. Next we have maintained that the so-designed spacetime is deducible 
from a modified non- Hilbert- Einstein field lagrangian, f(R) . Knowing the metric, 
we have solved the modified field equations retroactively for the sought-after f(R). 
Finally, we have compared our results with those of a conventional FRW model 
and have attributed the differences between the two to a dark energy component. 
Eventually, we have extracted the density, the pressure, and the equation of state 
of this stipulated energy. 

We note that our choice of the scale factor and the adjustment of its free parameters, 
to comply with the available cosmological observations, is, by no means, unique. 
The goal is simply to demonstrate that the use of the observations at the outset, 
to deduce the rudiments of what seems reasonable, facilitates the access to possible 
formal underlying theories, the action based f(R) gravity in our case. With the 
availability of more extensive and more accurate data in future one may come back 
and revise the model. See also for a similar emphasis. 

One of us (YS j^SI nas followed the same path to propose a modified gravitation for 
galactic environments and to explain the flat rotation curves and the Tully-Fisher 
relation in spiral galaxies without recourse to hypothetical dark matters. 
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Emprical a(x) 
from observation 




differential 
equation for F(x) 






boundary 
condition for F(x) 














R =R(X) 





Eliminating x 
results in f=f(R) 



equation of 
f=f(F, F',x) 



Fig. 6. Algorithm of inverse f(R): Choose a scale factor; Solve field equations first for F = df/dR 
and next for / as functions of the time parameter r; Eliminate r between R(t) /(t) to arrive at 

/(«)■ 
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